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CHAPTER I 


INTRODUCTION 


The problem of Fixed Satellite Service (FSS) system synthe- 
sis can be described as optimally allotting locations, fre- 
quencies and polarizations to communication satellites. 
These satellites are deployed in geostationary orbit (5.6 
earth radii above the equator) for transmitting signals to 
earth receivers. An orbital arc of feasible locations is 
determined for each satellite, based on the location and 
geometry of its service area. These arcs are comprised of 
locations from which every point in the intended service 
area is visible. A set of test points on the boundary of 
each sevice area is defined to facilitate the evaluation of 
a synthesis solution on the basis of interference. A fea- 
sible frequency band at which signals may be transmitted is 
specified, an accompanying polarization scheme is included 
as well. 

The primary objective in satellite synthesis is to allot 
orbital positions, frequencies and polarizations to satel- 
lites so that interference from unwanted satellite signals 
does not exceed a specified threshold at any (test) point 
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in any service area. The interference— geometry , based on 
single-entry interference between two down-link communica- 
tion circuits is shown in Fig. 1 111]. The following nota- 
tion is used: S - satellite, W - wanted network, I - inter- 
fering network, E - earth station, R - receiving antenna. 
In this figure, beams transmitted from the two satellites 
are aimed at their respective service areas. However, some 
signals transmitted from sateliite SI are received at the 
service area of satellite SU. These signals interfere with 
the signals transmitted from satellite SU and thus cause a 
deterioration in the quality of the signal of satellite SU 
received in its service area. The synthesis problem is to 
specify locations, frequencies and polarizations for satel- 
lites SU and SI so that interference between these two sat- 
ellites does not exceed a specified threshold anywhere in 
the service areas of satellites SU and SI. 

In real synthesis problems, there are typically many 
more than two satellites and two service areas. Aggregate 
c ar r i e r- to- in ter f e renc e (C/I) ratios calculated at every 
test point of each service area are used to measure the 
relative strength of wanted and interfering signals. Any 
solution in which all aggregate C/ I ratios exceed a certain 
threshold is considered acceptable. Uith this primary 
criterion in mind, other objectives for the FSS satellite 
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(a) Overall geometry 


sw SI 



(b) Detail with Earth radius exaggerated for clarity 


Figure 1 : INTERFERENCE GEOMETRY 
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system synthesis problem can be defined. As an example of 
this, we may wish to minimize the sum of the absolute devi- 
ations of allotted satellite locations from certain predet- 
ermined "desired" locations. 

The practical importance and political ramifications of 
a problem of this nature are obvious when one considers the 
proliferation of space technology to even the smallest and 
poorest countries. The need for orbital space and frequency 
and polarization allotments for present and future satel- 
lites illustrates the importance of prudent use of the 
geostationary orbit and the frequency spectrum. Since com- 
plex resource allocation problems are often attacked by 
formulating them as mathematical programming (optimization) 
problems, mathematical programming models and optimization 
methods described here should provide tools for the complex 
orbit/ spec trum allocation decisions to be made by those 
responsible for the governance of international communica- 
tions, such as delegates to the World Administrative Radio 
Conference's (WARC) and the International Telecommunication 
Union. 

In the present work, alternate mathematical programming 
formulations for FSS satellite synthesis are considered. 
Most of the models are similar but have different objective 
functions. Our aims are (1) to study different solution 
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techniques for satellite synthesis, and (2) to suggest an 
appropriate mathematical programming model, and an accompa- 
nying solution strategy, that can be reliably solved for 
synthesis solutions in a reasonable amount of time. 

The remainder of this manuscript is organized as fol- 
lows. Chapter II contains a review of literature on the 
subject of satellite synthesis. A nonlinear programming 
formulation of a satellite synthesis problem and three 
applicable solution methods are described in Chapter III. 
Nine mixed integer programming (MIP) formulations are 
described in the next chapter and shown in Appendix A. 
Results of computer runs with these nine models and five 
sma 1 1- to-med ium size geographically compatible scenarios 
are also presented and evaluated. In Chapter V, we present 
solutions to two of the models studied that were found with 
a heuristic solution procedure. Results for this approxi- 
mate method are compared with those obtained with branch- 
and-bound, an exact method. Solutions to a few large 
synthesis problems which have also been found with this 
heuristic are presented. Finally, in the last chapter, 
conclusions and recommendations about strategies for solv- 
ing satellite synthesis problems are discussed. 



CHAPTER II 


LITERATURE REVIEW 

MODELING APPROACHES FOR SATELLITE SYNTHESIS 

Satellite synthesis problems have been the subject of much 
research. Much of the work pertains to satellites in the 
Broadcasting Satellite Service (BSS). Solution strategies 
suggested for these problems range from exact algorithms to 
heuristic techniques. Some of these approaches can be mod- 
ified to suit the FSS also. 

There have been approaches that consider only frequency 
allotments. Cameron [3] has formulated an integer program- 
ming frequency allotment model that minimizes the number of 
channels allotted to satellites, subject to co-channel 
interference restrictions. Levis et al. [9] have suggest- 
ed another integer programming formulation for this same 
problem. They have also formulated an integer programming 
model that considers the allotment of multiple channels to 
each service area and takes into account adjacent-channel 
interference. Baybars [1] has formulated an integer 
programming model that minimizes the number of channels 
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used while considering both co-channel and adjacent-channel 
interference . 

Ito et al. [8] have formulated a satellite synthesis 
model that considers the allotment of satellite locations 
only. Their objective is to minimize the distance between 
the easternmost and westernmost allotted satellite loca- 
tions. The model is a nonlinear program that they suggest 
be solved using the Sequential Unconstrained Minimization 
Technique. In solving the model, satellites are positioned 
one at a time according to a predetermined launch sequence. 
Every satellite that is "launched" is positioned in the 
best possible manner relative to the other satellites that 
have already been launched. Such a method does not guaran- 
tee global optimality. 

Levis et al. [10] have formulated a nonlinear program- 
ming model for allotting locations and frequencies to BSS 
satellites. An alternating polarization scheme is assumed. 
Their model attempts to maximize the minimum of the aggre- 
gate C/ I ratios computed for all test points in every ser- 
vice area; locations and frequencies are restricted to lie 
within specified bounds. Levis et al. [10] and Martin et 
al. [14] have suggested that this model be solved using an 
extended gradient search procedure. Reilly et al. [18] 
have recommended a cyclic coordinate search procedure for 
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the same problem. An extensive computer experiment to 
evaluate the quality of the synthesis solutions found by 
these two search techniques was conducted by Reilly et al. 
[19]. They conclude that the cyclic coordinate search pro- 
cedure finds better solutions, but it does so at greater 
computing expense. This experiment and the two search 
methods are described in more detail in the next chapter. 

Heuristic methods for location, frequency and polariza- 
tion allotments have been suggested by Nedzela and Sidney 
[16], Chouinard and Vachon [4], and Christensen [5]* 


REQUIRED SATELLITE SEPARATION CONCEPTS 

Wang [22] has developed a procedure to calculate the 
required orbital separation between two satellites with a 
known mean longitude that guarantees single-entry co- 
channel Cl I ratios at test points along the boundaries of 
the intended service areas are at least equal to some 
acceptable threshold. Aggregate C/ I ratio requirements can 
be satisfied by specifying a higher threshold for the 
single-entry C/ I requirements [22]. 


We define is k to be the required separation between 
satellites i and j when their mean longitude is k degrees. 

Typically A<p values vary according to a bath-tub curve 
over the set of feasible locations [22]. 


We refer to the 
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maximum over all feasible k of the A <p . 


ijk 


as As.. 


The 


significance of As.^ is t ^ iat it is a worst-case required 
separation. Any pair of satellites i and j, when separated 


by As^, would satisfy our single-entry interference cri- 
teria irrespective of their locations in the feasible arc. 


In all the synthesis examples to be presented, we use 

^ 8 ij values calculated to ensure that aggregate C/ I 
requirements of 25dB are met. The A values are calcu- 
lated based on a single entry co-channel C/ I ratio require- 
ment of 30dB. The notion of requiring an extra 5dB over the 
single-entry interference criterion to meet an aggregate 
interference criterion was used at WARC '77 [6] and has 
proven to be valid for scenarios considered by Levis et al. 
[11] • The A0. values and the As*. values used in the 

1 J K 1 J 

experiments described in the forthcoming chapters have been 
calculated using Wang's [22] method. 


Detailed descriptions of the separation concept for sat- 
ellites with elliptical-beam antennas can be found in Wang 
[22], Levis et al. [11], and Reilly et al. [18]. Yamamura 
and Levis [21] have done similar work for satellites with 
circular-beam antennas, which is also described in Levis et 
al. [11]. 


Some of the models described in the literature have used 
similar separation concepts. A model developed by Ito et 
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al. [8] uses a separation matrix to ensure that single- 
entry interference requirements will be met. Required sat- 
ellite spacing calculations are based on the current loca- 
tions of the satellites already positioned ("launched"). 
Christensen [5] calculates separation values under the 
assumption that the required separation is approximately 
constant and independent of the satellite orbital longi- 
tudes. 

SYNTHESIS MODELS USING WANG * S SEPARATION CONCEPT 

The satellite separation concept has opened new avenues for 
the solution of satellite synthesis problems, as will be 
discussed in Chapters IV and V. The repeated, complex 
interference calculations seen in nonlinear programming 
models for satellite synthesis can be avoided. Such models 
have been used by Levis et al. 110], Martin et. al. 114], 
Reilly et. al. [19] and Ito et al [8]. Models that utilize 
the separation concept are a potentially valuable alterna- 
tive to the nonlinear programming synthesis models suggest- 
ed to date . 

Reilly et. al. [20] have incorporated the satellite sep- 
aration concept in a mixed integer programming (MIP) model 
and an almost linear programming model for allotting 
satellite locations. The logical relationships between 


satellite locations that are required to enforce the 
satellite separation constraints prohibit formulating the 
model as a pure linear program. Instead the simplex method 
with restricted basis entry, which does not ensure optimal- 
ity, is used to find solutions to the almost linear pro- 
gramming model 120]* The MIP model can be solved to opti- 
mality with a br anch-and-bound algorithm, but solution 
times for large problems may be prohibitive. The objective 
in both of these models is to minimize the sum of the devi- 
ations between prescribed (optimal) locations and the cor- 
responding "desired" satellite locations. Additional MIP 
formulations that use the satellite separation concept are 
described in Chapter IV. 

Gonsalvez [7J has developed a switching heuristic, or 
permutation algorithm, that is applicable to satellite 
synthesis problems as well as to some scheduling problems. 
He uses an MIP synthesis model suggested by Mo un t -Camp b e 1 1 
et al. [15] which incorporates the minimum required satel- 
lite separation concept. An initial ordering of the satel- 
lites is selected. With this initial ordering the synthe- 
sis problem reduces to a linear program. The linear program 
is then solved to find optimal satellite locations for the 
given ordering. The algorithm proceeds by considering per- 
mutations of small groups of adjacent satellites and re- 
optimizing the linear programming solution as needed. 
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Finite termination is guaranteed, but optimality is not. 
Feasible solutions to large synthesis problems (up to 81 
satellites) have been obtained with this heuristic in very 
reasonable solution times. Results obtained with this heu- 
ristic for some large problems are described in Chapter V. 
We also compare results for the switching heuristic with 


results for br anch-and-bound . 


CHAPTER III 


AN EVALUATION OF THREE SEARCH METHODS FOR A 
NONLINEAR PROGRAMMING FORMULATION FOR SATELLITE 

SYNTHESIS 

NONLINEAR PROGRAMMING FORMULATION 

As mentioned in the previous chapter, a nonlinear program- 
ming formulation of satellite synthesis problems has been 
suggested by Levis et al. [10]. The model as presented in 
Re illy et dl . [19] is: 

Minimize f(x) - I I I exp(a- (C/I)..,) (3.1) 

i jCJi kCKi ljk 

subject to, 

*i < x i i < w i Vi (3.2) 

l i ^ x f i i h i Vi (3.3) 

where x^ is the location of satellite i, x^ is the first 
1 i i 

(lowest) frequency in the family of channels assigned to 

satellite i, e ^ (w^) is the easternmost (westernmost) fea- 
sible location for satellite i, 1^ (h^) is the lowest 

(highest) frequency that can be assigned to satellite i, O 
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is a scaling constant, J ^ is the set of test points in the 
service area served by satellite i, and is the set of 
frequencies at which signals are to be transmitted from 
satellite i [19]. 


The objective function (3.1) is intended to maximize the 
worst aggregate C/ 1 ratio at any test point. By raising e 
to the (-C/lK^ power and summing over i,j, and k, the 
worst C/ 1 ratios contribute the most to the objective func- 
tion value. The constraints (3.2) and (3.3) are simple 
bound constraints and ensure that the location and frequen- 
cy decision variables remain feasible. 


This synthesis formulation is a difficult optimization 
problem to solve. The objective function in this model is 
not convex. In fact, it is characterized by numerous peaks 
and valleys of varying heights and depths [19J. The compu- 
tation of the C/ I ratios is much more complicated than it 
might appear from the formulation. The C/ I ratios are actu- 
ally a function of topocentric and satellite-centered 
angles, frequency discrimination, antenna gains and dis- 
crimination patterns, elliptical beam parameters, and 
transmitted power [19]. For brevity the entire expression 
for C/ I has not been presented. In this model C/ I is treat- 
ed as a function of locations and frequencies only. 
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COMPUTER EXPERIMENTS WITH THE GRADIENT AND CYCLI C 
COORDINATE SEARCH PROCEDURES 

Reilly et al. [19] have used the gradient search technique 
and the cyclic coordinate search technique to solve the 
model given in the previous section. 

The gradient search technique is a method that can be 
used to solve unconstrained nonlinear optimization prob- 
lems. It can be modified to solve problems with con- 
straints such as the ones present in the model we consider 
here. The method works as follows. We compute the gradient 
of the objective function at a feasible point. A line 
search is conducted in the negative gradient direction to 
find an improved solution* The method proceeds iteratively 
by restarting at the improved solution and computing the 
gradient again. The method terminates if a line search 
does not yield an improved solution, if a global optima 
(for convex programming problems) or a local optima (for 
non-convex programming problems) is found, or if the proce- 
dure goes through a predetermined number of iterations. 
Only solutions which are feasible are examined during the 
line search. 

The cyclic coordinate search procedure begins with the 
selection of a feasible point. Line searches are conduct- 
ed, in turn, in both the positive and negative coordinate 
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directions for each of the decision variables- The method 
terminates when a line search in the coordinate directions 
for each of the decision variables yields no improvement in 
the objective function value- The disadvantage of this 
method is that the objective function is evaluated many 
times. However, complex gradient computations are never 
done. 

Reilly et al. [19] designed an extensive computer exper- 
iment to assess the performance of the gradient and cyclic 
coordinate search algorithms as synthesis tools. A synthe- 
sis problem that consisted of seven South American coun- 
tries (Argentina, Bolivia, Brazil Chile, Paraguay, Peru 
and Uruguay) was used. Each country is to be served by one 
satellite. Each of the satellites is assumed to be capable 
of transmitting signals over a family of three channels 
with alternating polarization. Decision variables need 
only be defined for the lowest frequency assigned to each 
satellite. The remaining two frequencies are then the next 
two higher frequencies. All seven satellites are assumed 
to have the same feasible arc and feasible spectrum. 

Since the performance of the gradient and cyclic coordi- 
nate search techniques can be affected by various factors, 
a list of factors that could possibly affect performance 
was compiled. The factors and factor levels used are 



listed in Table 1. A 1 / 4 - f r ac t iona 1 factorial design, 

which included 64 runs, was chosen by Reilly et al. [19]. 

The important conclusions from the experiment are: 

(1) Thirty-two of the 34 best results were obtained by 
using the cyclic coordinate procedure. The average 
worst C/I for the gradient runs was 19.73dB, and for 
the cyclic coordinate runs the average was 41.10dB. 
Only three of the gradient search runs resulted in a 
worst C/I value greater than the acceptability thresh- 
old of 30dB. On the other hand, all cyclic coordinate 
runs satisfied this criterion of acceptability. 

(2) A cyclic coordinate iteration requires about 5.6 times 
as much CPU time as a gradient search iteration 
requires. (This conclusion is valid for the methods as 
they were implemented by Reilly et. al. [19]). 

(3) Better results are found when a longer feasible arc 
segment (80° - 110° W) is used. 

(4) When the starting locations are centered in the feasi- 

ble arc, the initial locations have little effect on 
the average of the observed worst C/ I ratios. However, 
when the satellites are located near the western bound- 
ary of the feasible arc, a better solution is found 
when the satellites are collocated than when they are 
separated by 1°. This is especially true for the 
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Table 1 

FACTORS AND FACTOR LEVELS 

Factor Levels 


Factor 

Low ( - 1 ) 

High (+1) 

A - Algorithm 

Gradient Search 

Cyclic Coordinate 
Search 

(or) 

Zoutendijk's Method. 

B - Location 
Sp ac ing 

0° 

1° 

C - Frequency 
Spac ing 

0 Mhz 

5 Mhz 

D - Starting 
Loca t ions 

Centered in 
Feasible Arc 

Spaced From Western- 
most Feasible Loc. 

E - Starting 
Frequency 

Centered in Ava- 
ilable Spectrum 

Spaced From Highest 
Available Spectrum 

F - Arc Length 

90° - 110° W 

80° - 110° W 

G - Frequency 
Spectrum 

12233-12300 Mhz 

12200 - 12300 Mhz 

H - Run Time 

5 CPU minutes or 
10 iterations or 
cycles 

10 CPU minutes or 
30 iterations or 
cycles 


Note. For the work carried out previously [19] a level of 
+1 for factor A represented the cyclic coordinate 
search procedure. For the present work a level of +1 
represents Zoutendijk's method. 
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gradient runs. 

(5) The poor performance of the gradient search algorithm 
can be attributed to a ,, jamming ,, effect that is often 
observed close to a boundary of the feasible region. 
This can be explained as follows. Some components of * 
the gradient may indicate that moves near a boundary 
look promising. The result will be that the line seg- 
ment searched in the negative gradient direction may be 
quite short, thereby prohibiting substantial reposi- 
tioning of the satellites. In effect, the algorithm 

becomes “jammed" against a boundary. 

COMPUTER EXPERIMENTS WITH ZOUTENDIJK'S METHOD 

The discussion above indicates that although better solu- 
tions are available, the "jamming" effect close to a bound- 
ary can result in poor performance of the gradient search 
algorithm. If we could avoid being trapped along a bound- 
ary, we could possibly find better solutions than those 
found with the gradient search with much shorter solution 
times than those obtained by the more computationally 
demanding cyclic coordinate search procedure. An alternate 
method to the gradient and cyclic coordinate searches, in 
terms of "jamming" and solution time, respectively, is 


Zoutendijk's method. 
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Zoutendijk's method, or the method of feasible direc- 
tions, is an algorithm which can prevent "jamming" [24]. 
The method involves solving a linear program to determine 
the most favorable direction in which to conduct a line 
search. The directions generated are a compromise between 
lining up with the negative gradient, as the gradient 
search does, and avoiding the boundary of the feasible 
region. For convex programming problems, this method will 
converge to a global optimum. However, as mentioned earlier 
in this chapter, the objective function in our problem is 
not convex and we would most likely find a local, rather 
than a global, optimum. At every iteration of Zoutendijk's 
algorithm a linear program has to be solved. Due to the 
special nature of our problem, this linear program is very 
sparse (Table 2) and reasonable solution times to the lin- 
ear program are usually observed. 

Given a feasible starting solution x to the nonlinear 

o 

optimization problem, 

Minimize z « f(x) (3.4) 

subjectto, 

g(x) < 0 (3.5) 

the d i r ec t ion- f ind i ng linear program to be solved is, 
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Table 2 

CONSTRAINT COEFFICIENT MATRIX OF SUBPROBLEM IN ZOUTENDIJK'S 

METHOD 


-V f +V f -V f +v 




+ 1 


-1 

+ 1 

0 

0 

0 

0 

0 . 


. . 0 

0 

+ 1 

+ 1 

-1 

0 

0 

0 

0 

0 . 


. . 0 

0 

+ 1 

0 

0 

-1 

+ 1 

0 

0 

0 . 


. . 0 

0 

+1 

0 

0 

+1 

-1 

0 

0 

0 . 


. . 0 

0 

+ 1 

0 

0 

0 

0 

-1 

+ 1 

0 . 


. . 0 

0 

+ 1 

0 

0 

0 

0 

+1 

-1 
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Minimize P (3.6) 

subject to, 

Vf( X q ) y - p <0 (3.7) 

g(x Q ) + Vg(x Q )y - pu < 0 (3.8) 

I |y.| < 1 (3.9) 

where p is a scalar, Vf(x Q ) is the gradient of the objec- 
tive function at x , Vg(x ) is the gradient of the constr- 

o o 

aints at x^, g(x^) is the value of the constraint functions 

(3.5) at x q , u — (1,1,1. ...,1), and y*(y 1> y 2 >y Q ) the 

direction vector. We replace | y ^ | with (y^ + + y^”)> where 
y^ + and y^ are restricted to non-negative values, to make 
constraint (3.9) linear before the linear program is solved. 

The algorithm for Zoutendijk's method is initiated at a 
feasible point. The linear program, given by equations 

(3.6) - (3.9), is solved to give us an optimal direction y 

in which to conduct a line search. Every time an improved 

solution is found, we start from the improved solution and 

solve another linear program. We have implemented the 

algorithm so that it terminates either when the solution to 

the linear program is zero, when 13 searches of 10 points 

each in the y direction yield no improved solution, or when 

m iterations are completed. The maximum number of itera- 
. * 

tions, m , was set at ten for the experiments we present 


here. 
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Formally the entire algorithm can be stated as follows; 


Step 


Choose a feasible starting solution x 

o 

Let in = 0. 



Step 2: If m ^ m , go to Step 6. Otherwise, compute 

f( ,X f ) . 

o 1 o 

Step 3: Solve the d irec t ion- f ind ing linear program. If the 

solution to this linear program is zero, go to 
Step 6. Otherwise, let n * 0 and m ■ m+ 1 . 


Step 4: Let n * n + 1 . Determine the length d of the line 

segment from (x^ ,x^ ) to the nearest boundary 

o o 

in the y direction. Evaluate the objective function 
(3.1) at ten equally spaced feasible points in the 
y direction, 

(*l,Xf) * (x x ,x f ) + kd . y , k * 1,2 10. 

o o 

(0.2) n 


Select the most favorable of these ten points, 


* * x 

( , x f ) . 

If f ( x ^* , x ) < f ( x ^ , 

o 


Xf ) 


/ \ / * * \ 
set ( x ^ , x £ ) = ( X£ ,x f ) 

1 o o 


and go to Step 2. Otherwise go to Step 5. 


Step 5: If n > 12, go to Step 6. Otherwise go to Step 4. 

Step 6: Stop. 
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The computer code used for the gradient runs in the 
experiments conducted by Reilly et al, [19] has been 
modified. FSS antenna patterns have been used instead of 
the BSS patterns used in [19]. All channels are assumed to 
have the same polarization. To solve the direction-finding 
linear program at each iteration, a linear programming 
package called ZX3LP, available through the IMSL Library, 
was linked to this computer code. For an n-satellite prob- 
lem, the linear program to be solved has 2n+l variables and 
4 n + 2 constraints - 4n constraints resulting from equation 
(3.8) above and one each for equations (3*7) and (3*9). The 
first constraint, resulting from (3.7) has coefficients 
equal to the gradient of the objective function with 
respect to each of the decision variables. The values of 
these gradient components range from 10~^ to 10~^. The 
rest of the constraint coefficients are 0 or 1. As one can 
see from Table 2, the constraint coefficient matrix of the 
linear program to be solved is very sparse. Due to the ill- 
conditioned constraint coefficient matrix, the gradient 
component values are scaled to avoid accuracy errors. 

A total of fourteen runs, seven each with the gradient 
search algorithm and Zoutendijk's algorithm, have been 
made. The same factors and factor levels used in [19] have 
been used here also. They are shown in Table 1. For the 
present work, a level of +1 for factor A (algorithm) 
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indicates Zoutendijk's method, and a level of -1 indicates 
the gradient search method. Factor H (run t ime/ i terat ion 
count) was set at the low level of -1 for all the runs. The 
results of this experiment are summarized in Table 3. A 
stop code of 1,2, or 3 indicates that the maximum iteration 
count was reached, that the allotted CPU time had expired, 
or that all attempts to find a better solution than the 
final reported solution failed. Finally, a stop code of 4 
indicates that the current solution is at a Kuhn-Tucker 
point, a local optima. (Note, however, that this did not 
occur) . 

All the runs summarized in Table 3 were made on an IBM 
3081 -D at The Ohio State University. 

CONCLUSIONS ON SEARCH METHODS FOR NONLINEAR PROGRAMMING 
FORMULATION 

It is seen from Table 3 that in each pair of runs, either 
Zoutendijk's method performed better than the gradient 
search or it yielded an acceptable solution, where accepta- 
ble means any solution with a worst C/I value of at least 
30dB. Only two of the gradient runs and three of the Zout- 
endijk runs meet our criterion of acceptability. In addi- 
tion to these, two other Zoutendijk runs (lb and 3b) came 
very close to the 30dB acceptability threshold. Both of 
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Table 3 

SUMMARY OF COMPUTER RESULTS 


Worst Iter./ CPU Stop 
C / 1 Cycles Time Type 


Run 

A 

B 

C 

D 

E 

F 

G 

H 



(sec) 


la 

-i 

+i 

+ i 

+ i 

-i 

+i 

-1 

-i 

5.59 

6 

264 

3 

lb 

+i 

+i 

+i 

+i 

-i 

+i 

-1 

-i 

29.03 

10 

77 

1 

2a 

-i 

+ i 

+i 

-i 

+i 

+i 

-1 

-i 

28.38 

. 5 

62 

3 

2b 

+i 

+i 

+ i 

-i 

+i 

+i 

-1 

-i 

36.05 

10 

64 

1 

3a 

-i 

+i 

+ i 

+ i 

+i 

-i 

+ 1 

-i 

5.59 

3 

122 

3 

3b 

+ i 

+i 

+ i 

+i 

+i 

-i 

+ 1 

-i 

24.39 

7 

85 

3 

4a 

-i 

+ i 

-i 

-i 

+i 

-i 

-1 

-i 

8.62 

2 

76 

3 

4b 

+ i 

+i 

-i 

-i 

+i 

-i 

-1 

-i 

17.21 

3 

88 

3 

5a 

-i 

+ i 

-i 

+i 

+i 

+i 

+ 1 

-i 

5.43 

8 

300 

2 

5b 

+ i 

+ i 

-i 

+ i 

+i 

+i 

+ 1 

-i 

29.41 

10 

61 

1 

6a 

-i 

-i 

+ i 

-i 

-i 

+i 

+ 1 

-i 

41.21 

5 

34 

3 

6 b 

+i 

-i 

+i 

-i 

-i 

+i 

+ 1 

-i 

32.62 

.10 

54 

1 

7a 

-i 

-i 

+i 

+ i 

-i 

-i 

+ 1 

-i 

47.97 

10 

32 

1 


4 7 


7b +1 -1 +1 +1 -1 -i +i -i 35 . 8 I 


10 


1 



these runs terminated because the maximum iteration count 


was exceeded. It is possible that another iteration would 
have increased the worst C/I value to an acceptable level. 
It is also observed that whenever the gradient search per- 
formed poorly, Zoutendijk's method found much improved 
solutions. On the other hand, for gradient runs which pro- 
duced acceptable solutions, Zoutendijk's method produced 
less attractive, but acceptable, solutions. The average 
G/I value for the gradient runs is 20.39dB while for the 
Zoutendijk runs it is 29*22dB. 

Runs made with Zoutendijk's algorithm always resulted in 
more iterations. Five of the seven Zoutendijk runs termi- 
nated because the maximum iteration count was exceeded. 
This implies that more improvement may have been possible 
in these runs. Only one of the gradient runs terminated 
due to the iteration count. No Zoutendijk run terminated 
at a Kuhn-Tucker point, again implying that more improve- 
ment was possible. 

Another significant observation that is made from the 
results presented here concerns solution times. The average 
run time per iteration for the gradient runs is 26.08 sec- 
onds while the corresponding value for the Zoutendijk runs 
is 10.25 seconds. The total average time for the gradient 
runs is 127 seconds, while for the Zoutendijk runs it is 68 
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seconds. Thus, on the whole, Zoutendijk's method appears to 
work faster. It can be seen from Run 5 that while the 
gradient search run had two iterations less than Zouten- 
dijk's algorithm, it took nearly five times as much time. 
Such a result can be explained by the fact that the gradi- 
ent search can become jammed against a boundary and perform 
fruitless line searches over very short line segments, 
thereby evaluating the objective function many times for 
naught. 

It is also observed that if the satellites are separated 
by 1° and are located at the western boundary, then the 
gradient search performs poorly. This is presumably because 
the easternmost satellites cause a blocking effect on the 
other satellites. For the same scenario, Zoutendijk's 
method, because of the nature of the algorithm, does very 
well. Another observation, that is not evident from the 
results presented here, is that satellite locations have a 
much greater effect than frequencies in determining the 
magnitude of the worst C/I values [2]. 

To summarize the findings from the experiments presented 
here it can be concluded that Zoutendijk's method looks 
very promising. It has performed better than the gradient 
search, both in terms of the worst C/I values and in terms 
of solution times. However, neither method guarantees an 
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acceptable aolution. To do full justice to the comparison 
of the gradient search and Zoutendijk algorithms, an 
extensive experiment, with several test problems, would 
have to be conducted. Alternate solution techniques, some 
of which are described in the next two chapters, are avail- 
able. These methods, unlike the nonlinear programming 
solution techniques, guarantee that any solution found is 
acceptable. 



CHAPTER IV 


MIXED INTEGER PROGRAMMING FORMULATIONS FOR 
SATELLITE SYNTHESIS PROBLEMS 

INTRODUCTION 


In this chapter we consider mixed integer programming (MIP) 
formulations of satellite synthesis problems. It is evi- 
dent from the discussion of the previous chapters that the 
major criterion to be satisfied in the satellite synthesis 
problem is that C/I ratios must not fall below a certain 
threshold value. The minimum satellite separation concept 
[22] is a valuable tool in ensuring that we do not exceed 
predetermined acceptable interference limits, and has been 
used in the formulations discussed below. This minimum 
satellite separation concept can be used in many synthesis 
models with a variety of objective functions. 

In the present work nine MIP models with different 
objective functions have been considered. The results of 
computer experiments with these models and five geographi- 
cally compatible scenarios is presented. Both point and 
arc allotments have been considered. The purpose of this 
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work is to evaluate the results from our different models 
in terms of robustness of objective function and solution 
quality. We try to identify objectives which are “univer- 
sal 11 in the sense that while meeting a particular criterion 
they do well against other objectives also. We look at the 
number of feasible solutions found with each model during a 
certain period of execution. Percentage improvements in 
objective function values are also considered. The initial 
ordering, the final ordering, and the change in the order- 
ing of satellites from the first feasible solution to the 
final solution are studied. 

OBJECTIVES FOR THE KIP FORKULAT IONS 


The nine different formulations (objectives) that have been 
cons idered are : 

1. Maximizing the minimum of the extra separation beyond 
A s, . for any pair of satellites i and j. This objective 
would tend to maximize the minimum C/I ratio and is thus 
similar to the objective used for the nonlinear program- 
ming formulation discussed in the previous chapter. 

2. Maximizing the minimum gap between adjacent satellites. 
The use of this objective would leave large gaps between 
satellites, facilitating introduction of new satellites 
between existing satellites at a future time. 



32 


3. Minimizing the total deviation of allotted locations 
from given "desired" locations, subject to feasible arc 
restrictions. The use of this objective can enable 
administrations to place satellites close to some pre- 
determined desired locations. In the experiments pre- 
sented here, the desired location used for each satel- 
lite is the center of its visible arc. 

4. Minimizing the maximum deviation of an allotted location 
from its given "desired" location, subject to feasible 
arc restrictions. This model attempts to do the same as 
model 3 except that we look at the maximum deviation 
rather than the total deviation of allotted locations 
from desired locations. 

5. Minimizing the total weighted deviation of allotted 
locations from given "desired" locations. This model 
fulfills the same objective as models 3 end 4 above. The 
weight used for each satellite is inversely proportional 
to the length of the satellite's feasible arc. This 
weighting scheme, by giving greater weight to adminis- 
trations with smaller feasible arcs, would attempt to 
ensure a feasible solution. 

6. Minimizing the maximum deviation of an allotted location 
from the closest boundary of the associated satellite's 
feasible arc. The solution to a synthesis problem with 
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this objective would leave a large gap in the center of 
the available arc enabling easy insertion of future sat- 
ellites. 

7. Minimizing the length of the arc between the easternmost 
and westernmost allotted locations. This objective 
attempts to place satellites as close together as possi- 
ble and thereby tends to conserve the geostationary 
orbit . 

8. Maximizing the minimum length arc segment allotted to 

any satellite. This objective differs from objectives 1 
through 7 in that it allots an arc segment to each sat- 
ellite rather than a point. This formulation gives an 

administration the flexibilty to place their satellite 
anywhere in the allotted arc segment. Arcs are permit- 
ted to overlap if the minimum required satellite separa- 
tion is zero. This model allots each administration the 
same length arc segment. 

9. Maximizing the minimum length weighted arc segment 
allotted to any satellite. This objective might be more 
practical than that in model 8 above because we could 
allot arcs to administrations (countries) such that the 
length of each allotted arc segment is proportional to 
some criterion and also give administrations the 
flexibility to place satellites anywhere in their 
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allotted arc segments. The weighting criterion to be 
used could be a country's population) telephone traffic, 
land area, population density, etc. The weights we use 
here are proportional to the population of each country. 

For the sake of brevity, only the formulation for objective 
1 is given below* Formulations (2) through (9) and the 
notation used in them is given in Appendix A. The notation 
that we use for formulation (model) 1 is given below; 


Par ame ter s : 


vv 


easternmost (westernmost) feasible location for 
satellite j, in degrees west longitude. 


E * minimum l E j] over all j (assumed to be zero). 


W 


maximum [W.] over all j. 


As.. = worst-case minimum required satellite separation 
between satellites i and j, in degrees longitude. 


Decision variables: 


x . 
J 


orbital location of satellite j (degrees west of E), 
in degrees west longitude. 


x . . 
1 J 



if satellite i is located west of satellite j, 
o the rwi s e 


a = minimimum of the extra separation beyond As., for all 
pairs of satellites i and j. 
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Model 1 : Maximize the minimum of the extra separation beyond 
a required minimum satellite separation. 

Maximize z * a (4*1) 


subject to, 


( x.-Xj) + [2(E-W>- AsjjJxjj - a > 2(E-W) (4.2) Vi, j 3 i<j 

(-x i + X j) - [ 2 ( E- W ) - As-^Jx^ - a > As.. (4.3) ¥i, j 3i<j 

E ; < x j < w j (4.4) Vj 

x^ = 0 or 1 (4.5) V i , j 3i<j 


Xj > 0 (4.6) Vj 

a > 0 (4.7) 

The objective function (4.1) maximizes the minimum of the 
extra separation between two satellites beyond their 
required minimum separation. For each pair of satellites i 
and j, exactly one of the constraints (4.2) and (4.3) will 
be redundant. The binary variable x^ which takes value 1 
or o , depending on whether satellite i is located to the 
west of satellite j, determines the redundant constraint. 
If x £j“l then constraint (4.3) is redundant. If x^j«0 then 
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constraint (4-2) is redundant. The non-redundant con- 
straint then enforces the required separation between sat- 
ellites i and j. Constraints (4*2) and (4*3) also measure 
the minimum separation beyond a required satellite separa- 
tion, which is denoted by a. Constraints (4*2) and (4.3) 
are an example of dichotomous constraints commonly used in 
mixed integer programming formulations to enforce "either- 
or M conditions. 

Constraints of type ( 4 . 4 ) ensure that all allotted sat- 
ellite locations are feasible. In the M1P code used to 
solve this synthesis problem, constraint (4.4) is enforced 
as a simple bound on the variable x^. The integrality of 
the binary variables is enforced by (4*5) end non- 
negativity requirements of continuous variables is ensured 
by (4.6) and (4.7). 

For a problem with m satellites, this formulation 
entails m(ni-l) structural constraints (4.2) and (4.3), m + 1 
continuous variables and m(m-l)/2 binary variables. 

The formulations given in Appendix A are very similar to 
the one shown above. Each of these formulations uses the 
minimum required satellite separation, As^ . ( However, in 
several of the formulations, additional variables are 


introduced 



COMPUTER EXPERIMENTS WITH MIP MODELS 


Five scenarios were used to evaluate the nine MIP models. 
These scenarios are listed below. The countries and satel- 
lites which make up each scenario are given in Table 4. 

( 1 ) E. Europe (12 s a t e 1 1 i t e s / 1 2 service areas). 

(2) W. Europe (12 sate Hi tes/ 1 2 service areas). 

(3) S* America (13 s a t e 1 1 i t e s / 1 2 service areas). 

(4) N. Africa (10 s ate 1 li tes/ 1 0 service areas). 

(5) S • E . Asia (10 s a t e 1 1 i t es/ 1 0 service areas). 


The results of computer experiments with these five sce- 
narios and our nine models are given in Appendix B. An IBM 
linear p rog rammi ng /mixed integer programming package called 
MIP/370 has been used for all the runs. The runs for these 
experiments were made on an IBM 4341 at The Ohio State Uni- 
versity. The maximum CPU time specified was 20 minutes and 
maximum node storage was set at 5000 nodes. Computer runs 
were automatically terminated if either of these limits was 
exceeded . 

To evaluate the universality, or robustness, of each of 
the models, the final solution from formulations 1 through 
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Table 4 

GEOGRAPHIC SCENARIOS USED IN EXPERIMENTS 



EAST 

EUROPE 

WEST 

EUROPE 

SOUTH 

AMERICA 

NORTH 

AFRICA 

1 . 

Fi nland 

Italy 

Bolivia 

Libya 

2 . 

Bulgaria 

Norway 

Brazil* 

Nigeria 

3 . 

Romania 

Denmark 

Chile 

Mali 

4 • 

Greece 

Belgium 

Co lomb ia 

Morocco 

5 . 

Albania 

France 

Guyana 

Sudan 

6. 

Poland 

Switzerland 

Peru 

Egypt 

7 • 

Hunga ry 

Nether land 

Paraguay 

Chad 

8 • 

Swede n 

Spain 

Equador 

Tunisia 

9 . 

Austria 

Ire land 

Venezuala 

Algeria 

10. 

E. Germany 

W. Germany 

Argent ina 

Mauritius 

1 1 . 

Yugos lavia 

Por tuga 1 

Uruguay 


12 . 

Czechoz- 
1 avakia 

United 

Kingdom 

Surinam/ 

F. Guyana 



S. E. 
ASIA 


Taiwan 

Cambodia 

Malaysia 

Vietnam 

China 

Burma 

Thai land 

Laos 

Indones ia 

Philipp- 

ines 


* _ 


2 Satellites 


39 


7 was evaluated in the objective function of each of the 
other six models for each scenario. The results of this 
analysis can be found in Appendix C . The initial ordering, 
the final ordering, and the change in the ordering of sat- 
ellites from the first feasible solution to the final solu- 
tion was also investigated. Our findings are given in 
Appendix D. 

Table 5 gives the average number of feasible solutions 
found with each of the models for the five scenarios. The 
average percentage improvement and the average percentage 
improvement with respect to time for each of the models is 
also given in Table 5« 

From the results given in Appendices B,C and D, and in 
Table 5* we make the observations discussed in the next 
section. 

2®servations_of_cohputer_runs_on _the_mip_models 

As mentioned earlier in this chapter, we have considered 
both point allotment models (models 1 through 7) and arc 
allotment models (models 8 and 9)* We will analyze the 
results of computer runs of the point allotment models sep- 
arately from the results for the arc allotment models. 
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Point _Allotnent Models 

On an examination of the results given in the appendices 
and from the discussion below it will become clear that 
models 3,4? and 5 are robust. Furthermore, solution times 
for these models are less than those for the other point 
allotment models. Hence, in addition to overall results, 
we will also give results pertaining to these three models. 

The tables in Appendix B give us the times to the first 
feasible solution and the final feasible solution for each 
of the branch-and-bound runs. We see that for all our mod- 
els and for each of the scenarios considered, the first 
feasible solution was found very fast. For all the point 
allotment models, the average time to the first feasible 
solution was 11.9 seconds. On the other hand, the time to 
the first feasible solution for models 3,4, and 5 was only 
5.5 seconds. 

As mentioned earlier, the maximum CPO time specified for 
all the runs was 1200 seconds. However, we find little or 
no improvement in the objective function values during the 
last few minutes of execution, even though optimal solu- 
tions were found for but a few of the runs. The average 
time to the best solution for all point allotment models 
was 385 seconds. On computing the same quantity for models 
3,4, and 5 we get a value of 165 seconds. An optimal 


41 


solution was found and proved for models 3,4, and 5 in 
three of the five scenarios. A proven optimal solution was 
found for model 6 in one scenario. All other point allot- 
ment runs terminated at the maximum allotted CPU time. 

Table 5 gives us the average number of feasible solu- 
tions found for each of our models. It is seen that model 
2 found an unusually large number of feasible solutions. 
This can be explained as follows. Recall that model 2 max- 
imizes the minimum gap between adjacent satellites. Once 
the entire arc is occupied, the gap between each pair of 
adjacent satellites is the same. In each of our runs, this 
gap was larger than the largest required satellite separa- 
tion value. Hence, alternate solutions to the model are 
found by switching satellites from one position to another. 
The relative locations and hence the objective function 
values remain the same for such solutions. Each time a 
satellite is switched we find a different feasible solu- 
tion. Since such switches occur many times we find many 
different feasible solutions. 

We also see from Table 5 that models 3,4, and 5 find a 
relatively small number of feasible solutions. This is 
because some of the runs terminated when an optimal solu- 
tion was found, before the maximum CPU time expired. The 
numbers of feasible solutions found in 20 minutes of CPU 
time for models 1, 6 and 7, are nearly the same. 
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Table 5 

AVERAGE NUMBER OF FE ASIBLE SOLUTIONSAND PERCENTAGE 
IMPROVEMENTS IN OBJECTIVE FUNCTION VALUE 


AVG. NO. OF % IMPROVEMENT % IMPROVEMENT 

FEASIBLE PER CPU MINUTE 

SOLUTIONS 


Model 

1 

13.0 

6 4.16 

13 . 69 

Model 

2 

68.0 

31.48 

15.84 

Model 

3 

6.4 

20.20 

18.70 

Model 

4 

3.6 

11.35 

29.69 

Model 

5 

6.4 

25.41 

53.59 

Model 

6 

10.4 

82.38 

35.16 

Model 

7 

11.8 

21.32 

1 . 93 

Model 

8 

10.0 

38.64 

10.88 

Model 

9 

8.0 

22 . 80 

9.41 


NOTE. All the above runs were given a maximum CPU time of 
20 minutes. Some of the runs reached optimality 
before the maximum time expired. 
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Table 5 also gives the percentage improvement in solu- 
tion value from the first feasible solution to the best 
solution, and the percentage improvement with respect to 
the time to find the best solution, A significant point 
emerging from this analysis concerns models 3 and 5* 
Recall that model 3 minimizes the total deviation from 
desired locations. In model 5 also we minimize the sum of 
deviations but with a weighted objective function and with- 
out enforcing the visible arc restrictions. From Table 5 
we see that model 5 converges to its final solution faster 
than model 3* Appendix C, in which we compare the solution 
for each model against the objective function values of the 
other models, shows us that the solution for both of these 
models is similar. We thus conclude that it might be 
advantageous to use model 5 instead of model 3. However, 
we must note that the solution from model 5 might not be 
feasible, in which case we would have to use model 3. 

We see from Table 5 that models 1 and 6 show a signifi- 
cant improvement in objective function value from the first 
feasible solution to the best solution found. Model 4 shows 
the least percentage improvement from the first feasible 
solution to the best solution. We also note from Table 5 
that model 5 has the fastest convergence to the final solu- 
tion, while model 7 converges very slowly to its final 


solution. 
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The solutions from models 1 through 7 were evaluated in 
the other six models by calculating their "would-be" objec- 
tive function values* These solution values are reported 
in Appendix C. The following observations have been made: 

(1) The best solution to model 1 gave an objective function' 
value close to that of the best solution to model 2, 
and vice versa • 

(2) Also, the best solution to model 3 gave an objective 
function value close to that of the best solution to 
model 5, and vice versa. 

(3) The solutions to models 3,4, and 5 gave good values for 
the objective function of model 7, but the reverse was 
not true. Recall that in our case the desired locations 
are assumed to be the center of the feasible arcs. By 
minimizing deviations of prescribed satellite locations 
from desired locations, we also conserve the geosta- 
tionary orbit, but not vice versa. 

In models 8 and 9 we allot arc segments to satellites. In 
model 8 the length of the arc segment allotted to every 
satellite (administration) is the same. In model 9 the the 
length of the arc segment allotted to every administration 
is proportional to the population of the administration. 
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From Appendix B we see that the overall average time to the 
first feasible solution for model 8 is 11.0 seconds while 
for model 9 it is 7*5 seconds. The overall average time to 
the best solution found is 557 seconds for model 8 and 333 
seconds for model 9. 

A proven optimal solution was found for model 8 in one 
scenario and for model 9 in two scenarios. From Table 5 we 
see that the average number of feasible solutions found for 
models 8 and 9 is ten and eight, respectively. 

From Table 5 we also note that model 8 shows a greater 
percentage improvement in objective function value from the 
first feasible solution to the final solution than model 9. 
However, the rate of convergence for models 8 and 9 is 
nearly the same. 

Ana lysis of Satellite Orderings 

Appendix D gives the ordering of satellites for each model 
and each scenario. In this Appendix, ( F) refers to the 
ordering at the first feasible solution and (B) is the 
ordering of satellites at the best solution found. For the 
arc allotment models we use the center of the allotted arc 
to determine the order of satellite allotments. For all 
the point and arc allotment models considered we see a sig- 
nificant change in the orderings from the first feasible 
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solution to the best solution. It is also seen that, for 
the final solution found, there is some similarity in the 
ordering at the boundaries. For example, from Table 21 we 
see that satellites 4, 5, and 11 are found close to or at 
the boundaries of the arc occupied by the eastern European 
satellites. From Table 24 we see that satellites 5,6, and 
7 are located close to the boundary of the arc occupied by. 
the African satellites. Similar results are found for oth- 
er scenarios. Satellites which show this kind of similari- 
ty at the boundaries usually have feasible arcs whose east- 
ern and/or western edges extend beyond those for other 
satellites. Such observations may be useful in developing 
heuristic techniques which reduce the number of binary 
variables. Satellites with feasible arcs extending beyond 
the eastern or western edges of the arcs of other satel- 
lites could be fixed in position. As an example, for a 
13-satellite problem, if we could fix four satellites in 
positions at either end of the ordering, we would reduce 
the number of binary variables from 78 to 36. This signifi- 
cant reduction in the number of binary variables reduces 
the complexity of the problem considerably. 


CHAPTER V 


APPLYING A HEURISTIC ALGORITHM FOR SATELLITE 

SYNTHESIS 


The satellite synthesis problem can be considered to be 
comprised of essentially two problems - a satellite order- 
ing problem and a satellite location problem (Hount- 
Campbell et al. [15]) • For a given ordering of satellites, 
the synthesis problem reduces to a linear program. This 
linear program can be solved to find the optimal locations 
of satellites subject to required minimum satellite separa- 
tion and feasibile arc constraints. The determination of 
an optimal, or even good, ordering is still difficult. 

Gonsalvez [7] has developed a heuristic algorithm that 
uses this notion to find solutions to satellite synthesis 
problems. In the current implementation of this algorithm, 
we have the option to minimize either the sum of deviations 
of allotted locations from desired locations or the maximum 
deviation of an allotted location from the corresponding 
desired location [7j. 
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The heuristic algorithm is initiated with the specifica- 
tion of an initial ordering of satellites. Given this 
ordering we solve a linear program to find optimal satel- 
lite locations. Next, we generate another ordering by per- 
muting ( M s witching 1 *) the order of satellites within a group 
of k adjacent satellites. For every promising switch, the 
linear programming solution is reoptimized using the prin- 
ciples of sensitivity analysis. If a proposed switch does 
not yield an improved solution, the switch is reversed. If 
an improved solution is found, that solution becomes the 
incumbent and we look for other promising switches starting 
from this improved solution. This procedure is repeated 
till all orderings that can be obtained by permuting k 
adjacent satellites at a time have been considered. When 
no promising permutations are identified, k is incremented 
by 1 (to a maximum of 5), or the algorithm is terminated. 

To guarantee optimality of the solution found with the 
H swi tching" heuristic, one would have to examine all n! 
orderings, where n is the number of satellites. This is 
not practical for large problems. Since the objective is 
to minimize the sum of deviations from desired locations or 
the maximum deviation from a desired location, we begin 
with the initial ordering based on the order of desired 
locations and permute small groups of adjacent satellites 
to generate different orderings. 


4 9 


Formally the algorithm for the switching heuristic can 
be stated as follows; 

Step 1; Generate an ordering of the satellites based on 
their desired locations* 

Step 2: Given the ordering from Step 1, and a particular 

objective, solve a linear program to find the opti- 
mal locations for the satellites. 

Step 3: Select a particular ad j ac e n t - s a t e 1 1 i t e group size, 

(2<k<5) . 

Permute the order of satellites within this group 
(i.e., switch satellites in position 1 to k with 
the satellites in positions k + 1 to n being fixed in 
position (not location)). 

Reoptimize the linear programming solution as need- 
ed. 

Repeat till all orderings within a group have been 
considered . 

Step 4: Repeat Step 3, switching satellites 2 to k+1, 3 to 

k + 2 , n-k+1 to n. 

Step 5: (This step is optional). 

Increase group size by 1. Start from Step 3 with 
the best solution found so far. 

Step 6: If no improved solution is found, stop. 
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To solve the linear program, the algorithm uses the for- 
mulation of model 3 suggested by Mo un t -Camp be 1 1 et al. 
[15]. It has been found that increasing the group size k 
from 2 to 5, in steps of 1, is most effective [ 7 J • For 
further details on the switching heuristic the reader is 
referred to [ 7 ] • 

The switching algorithm can use either the worst-case 
required satellite separations (As) or the location- 
dependent satellite separations ( A</>) to guarantee inter- 
ference requirements are met. If the A0 values are to be 
used, then A <£ matrices, calculated for selected satellite 
longitudes, must be input. This option cannot be imple- 
mented conveniently in a b r anch-and-bound algorithm. It 
represents a distinct advantage of the switching algorithm, 
because the smaller location-dependent satellite separa- 
tions may make it easier to find feasible solutions to 
t igh t ly - c o ns t r a i ned problems. 

We have used the switching heuristic to find solutions 
to some large and medium-size problems, in addition to the 
five scenarios introduced in Chapter IV. In Table 6 
results for this algorithm on a 36-satellite scenario are 
given. Results obtained from a larger 59-satellite scenar- 
io are given in Table 7. Results from a 81-satellite scen- 
ario are given in Table 8. 
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In Table 9 we present results obtained with the switch- 
ing heuristic for the five scenarios used to evaluate the 
MIP formulations in the previous chapter. The objective 
function is to minimize the sum of deviations from desired 
locations. We compare these results with the results 
obtained for the MIP runs. The CPU time for the switching 
runs is the total run time. For the MIP runs, it is the 
time to the best solution if an optimal solution was not 
found or the total run time if optimality was proved. We 
must note that the CPU times are not exactly comparable 
since two different computer systems have been used. We do 
know that the IBM 3081-D is faster than the IBM 4341* 

Table 10 gives results similar to those in Table 9, but 
the objective function is to minimize the maximum deviation 
from given desired locations. 

From Tables 6,7, and 8 we see that the switching heuris- 
tic finds feasible solutions very fast. A 36-satellite 
scenario, comprised of European and African satellites, was 
modelled using formulation 3, the minimization of the total 
deviation from desired locations, given in Appendix A. We 
ran this model on the IBM 4341 using MIP/370 and found no 
feasible solution in 120 minutes of CPU time. On the other 
hand, the switching heuristic found a feasible solution to 
this same scenario in just 0«4 seconds. The times obtained 
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Table 6 

SWITCHING MMMlIIiLWITH 3 6 _S AT E L L I T E S 

OBJECTIVE: MINIMIZE THE TOTAL DEVIATION FROM GIVEN DESIRED 
LOCATIONS . 

CODE USED: SWITCHING HEURISTIC WITH INCREASING 'k'. 

SCENARIO : EUROPE + AFRICA (36 SATELLITES). 

DELTA_^S_^ 

“tot~dev7 max“dev7 orbital cpu“time 

ARC (Sec) 


FIRST FEASIBLE 642.49 38.17 100.45 0.40 

SOLN. 


FINAL SOLN. 169.29 25*32 72.56 436 


DELTA '<#>' 


TOT DEV. MAX DEV . ORBITAL CPU TIME 

ARC (Sec) 


FIRST FEASIBLE 642.19 38.17 100.45 1.60 

SOLN. 


FINAL SOLN. 121.79 14.08 72.00 541 
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Table 7 

SWITCHING HEURISTIC WITH 59 SA TE LLITES 

OBJECTIVE: MINIMIZE THE TOTAL DEVIATION FROM GIVEN DESIRED 
LOCATIONS . 

CODE USED: SWITCHING HEURISTIC WITH INCREASING 'k'. 

SCENARIO : OASTS2G1 (59 SATELLITES). 

delta 

TOT DEV. MAX DEV. ORBITAL CPU TIME 

ARC (Sec) 


FIRST FEASIBLE - 1200* 

SOLN. 


DELTA '0' 


TOT DEV. MAX DEV. ORBITAL CPU TIME 

ARC ( Sec ) 


FIRST FEASIBLE 475.55 26.95 94.15 153-8 

SOLN. 


FINAL SOLN. 443*98 24.54 92.35 1200* 


* - TERMINATING DUE TO TIME LIMIT 


54 


Table 8 

SWITCHING HE U RI STIC WITH 81 SA TELL ITES 

OBJECTIVE: MINIMIZE THE TOTAL DEVIATION FROM GIVEN DESIRED 
LOCATIONS. 

CODE USED: SWITCHING HEURISTIC WITH INCREASING 'k'. 

SCENARIO : OASTS2G1 + CARRIBEAN (81 SATELLITES). 





DELTA 

111 





TOT DEV. 

MAX 

DEV. 

ORBITAL 

ARC 

CPU TIME 
( Sec ) 

FIRST 
S OL N • 

FEASIBLE 




- 

600* 




DELTA 

'0' 





TOT DEV. 

MAX 

DEV. 

ORBITAL 

ARC 

CPU TIME 
( Sec ) 

FIRST 

SOLN. 

FEASIBLE 

1033.29 

53 - 

5 3 

133.07 

354 . 7 

FINAL 

SOLN. 

832.46 

47 . 

01 

126.55 

* 

600 

- 

— 










* _ 


TERMINATING DUE TO TIME LIMIT 
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for the 59-satellite and 81-satellite scenarios are also 
r easonab le . 

From Tables 9 and 10 we see that the switching heuristic 
does better than the MIP formulations in terras of solution 
time in nearly all cases. The objective function values 
are also comparable. in six of the ten runs given in these 
two tables, we find that the switching heuristic came up 
with the same solution as the MIP solution. In four of 
these cases the switching heuristic found an optimal solu- 
tion. In the other four computer runs, where we did not 
find the same solution with both approaches, the switching 
heuristic solution's value is nearly the same as that for 
the best solution found by the branch-and-bound algorithm. 
The results presented here indicate that the switching heu- 
ristic performs exceedingly well, both in terras of solution 
time and solution quality. The advantages of this algor- 
ithm can be summarized as follows: 

(1) The switching heuristic finds feasible solutions very 
fast. 

(2) The final solution values are nearly as good as those 
obtained by truncated runs with a branch-and-bound 
algorithm. In some cases, the switching heuristic con- 
verges to the optimum solution. 
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Table 9 

SWITCHING H EURISTIC VS. MIP MODE L 3 

OBJECTIVE: MINIMIZE THE TOTAL DEVIATION FROM GIVEN DESIRED 
LOCATIONS . 


CODE USED: SWITCHING HEURISTIC WITH INCREASING 'k'. 


OPTION ; DELTA 

'S' 





SWITCHING HEURe 

MIP FORMULATION 

SCENARIO 

Ob j . Fn . 

CPU Time 

Ob j . Fn . 

CPU Time 


Value 

(Sec) 

Value 

(Sec) 

E. EUROPE (12) 

52 . 74 

7.34 

49.87 

613.0 

W. EUROPE (12) 

32.29 

11.34 

★ 

29.67 

41 .4 

S . AMERICA (13) 

30.44 

11.35 

30.44 

86.4 

N. AFRICA (10) 

8.65 

3.68 

8.65* 

3.6 

S.E. ASIA (10) 

23 .05 

5-31 

23.05* 

96.0 


NOTE: SW runs have been made on the IBM 3081-D. 

MIP runs have been made on the IBM 4341 . 

For SW runs, CPU time = total run time 

For MIP runs CPU time = total run time for optimal 

solutions, (or) 

= time to best solution 
if soln. is not optimal 

Figures in parentheses indicate number of satellites. 


* _ 


PROVEN OPTIMUM 
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Table 10 

SW1TCHING_HEURISTIC_VS^_MIP_M0DEL_4 

OBJECTIVE: MINIMIZE THE MAXIMUM DEVIATION FROM GIVEN 
DESIRED LOCATIONS. 

CODE USED: SWITCHING HEURISTIC WITH INCREASING 'k.'. 


OPTION : 

DELTA 

'S' 






SWITCHING HEUR. 

MIP FORMULATION 

SCENARIO 


Ob j • Fn . 
Value 

CPU Time 
. (Sec) 

Ob j . Fn . 
Value 

CPU Time 
( Sec) 

E. EUROPE 

(12) 

9.41 

8.68 

8.48 

1137.0 

W. EUROPE 

(12) 

6.62 

8.35 

* 

6.62 

155.4 

S. AMERICA (13) 

6.96 

10.23 

6.60 

12.6 

N. AFRICA 

(10) 

1 .64 

3.29 

1 • 64* 

0.4 

S.E. ASIA 

(10) 

5.19 

5.10 

5.19* 

79 . 8 


NOTE: SW runs have been made on the IBM 3081-D. 

MI P runs have been made on the IBM 4341. 

For SW Runs CPU time ■ total run time 

For MIP runs CPU time ** total run time for optimal 

solutions 

* time to best solution 
if soln. is not optimal 

Figures in parentheses indicate number of satellites 


* _ 


PROVEN OPTIMUM 
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(3) The switching heuristic is independent of the objective 

functions The current implementation includes the 

option of choosing between two objective functions; the 
algorithm can be modified to solve other models also, 
including arc allotments [ 7 ] • 

(4) The switching heuristic can make use of location- 

dependent satellite separations (A <f> ) values. This 

means that it is more likely to find a feasible solu- 
tions when one exists, especially for tightly- 


constrained problems 



CHAPTER VI 


RECOMMENDATIONS 

Our aims were (1) to analyze different solution techniques 
for solving satellite synthesis problems and (2) to suggest 
an appropriate mathematical programming model, and an 
accompanying solution strategy, that can reliably solve 
synthesis problems at reasonable cost. In the previous 
chapters we have discussed three distinct solution strat- 
egies, namely, nonlinear programming search methods, 
branch-and-bound, and a switching heuristic. We have also 
experimented with different satellite synthesis models. We 
have evaluated these strategies and the models. 

In Chapter III we see that though it does not guarantee 
that C/I needs will be met, Zoutendijk's method does better 
than the gradient search in terms of solution time and 
solution quality. The "jamming 11 effect, which is a poten- 
tial problem with the gradient search method, is not 
observed when Zoutendijk's algorithm is used. Results pre- 
sented by Reilly et al. 119] indicate that the cyclic 
coordinate search method is likely to give worst C/I values 
that are better than those obtained with Zoutendijk's 
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method. However, solution times are likely to be much 
longer. It must be pointed out that even the cyclic coor- 
dinate search method, like the gradient search and Zouten- 
dijk's method, does not ensure an acceptable solution. 
Approaches that guarantee that acceptable solutions will be 
found, provided such solutions exist, are more useful. 

In Chapter IV we see that by introducing the minimum 
satellite separation concept into an MIP model we guarantee 
that any solution found will meet our interference require- 
ments and hence be acceptable. We have also seen that for 
the problems ve have considered here (10-13 satellites) the 
first feasible solution is found very fast. However, any 
real-life problem would be larger than this. Solution 
times for integer programming problems typically increase 
exponentially with problem size. An example of this was 
mentioned in Chapter V where we said that no solution was 
found to a 36-satellite MIP model in two hours of CPU time. 

The nonlinear programming methods discussed in Chapter 
III, unlike the branch-and-bound algorithm, do not ensure 
that globally optimum solutions will be found. We have 
seen that in spite of extensive line searches, none of the 
nonlinear programming algorithms converged to even a local 
optimum. On the other hand, a few of the MIP models have 
homed in on globally optimum solutions. For the smaller 
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( 7 -s a t e l 1 i t e ) problems solved by the nonlinear programming 
methods, the time to find the first feasible solution, if 
one is found, is greater than the time to find the first 
feasible solution for the larger (10 to 13-satellite) MIP 
models (2]. Also, the MIP models guarantee an acceptable 
solution if a feasible solution is found. Clearly, the 
mixed integer programming models are preferred over the 
nonlinear programming model of Chapter XI. 

The switching heuristic finds feasible solutions rela- 
tively quickly* Solution quality is also nearly as good as 
that obtained by truncated branch-and-bound runs. In addi- 
tion to this, the switching heuristic makes use of the 
smaller, location-dependent satellite separations ( A<f>) * 
It is independent of the objective function and can be mod- 
ified to solve models other than the deviation models we 
considered in Chapter V. 

From the experiments presented in the previous two chap- 
ters, it is clear that for small problems the solutions to 
models solved by the truncated branch-and-bound method are 
as good or slightly better than the solutions to similar 
models solved by the switching algorithm. If computer time 
is limited, then we recommend that the switching heuristic 
be used. For larger problems, the MIP models would require 
prohibitive amounts of computing resources. Such problems 
are better solved by the switching heuristic. 
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Since the nonlinear programming methods and branch-and- 
bound have known computational shortcomings, we conclude 
that the heuristic procedure is the best of the methods 
considered for solving satellite synthesis problems. 

Our other important conclusions come from Chapter IV in' 
which we compare mixed integer programming models with dif- 
ferent objective functions. In that chapter we show that 
the objective of minimizing deviations from desired loca- 
tions (models 3,4, and 5) have an edge over the other 
objectives because: 

(1) Optimal solutions are found with models 3,4, and 5 a 
greater number of times than with any other model. 

(2) The time to the first feasible solution and the time to 

the best solution for models 3,4, and 5 are less than 

half of the overall average times for the other point 
allotment models. 

(3) The solutions to models 3,4, and 5 give good values in 

the objective function of model 7 (minimizing the total 

arc occupied by the satellites to be allotted loca- 
tions) . 

Summarizing, we recommend that for point allotment 
problems model 3,4, or 5 be used. Model 5 might not find a 
feasible solution for some problems and one must keep this 
in mind while deciding between models 3 and 5. 
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The two arc allotment models we consider are models 8 
and 9. We have mentioned in Chapter IV that model 8 finds 
solutions in which the length of the arc segment allotted 
to every administration is the same. On the other hand, 
model 9 finds solutions in which the length of the arc seg- 
ment allotted to every administration is proportional to 
some criterion. We also note that the time to the first 
feasible solution and the time to the best solution is less 
for model 9 than for model 8. Hence, for arc allotment 
problems, model 9 seems more practical. 


MIXED INTEGER PROGRAMMING MODELS 
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NOTATION 


We define the following notation for the formulations dis- 
cussed in this Appendix: 


Parameters . 


E.(W.) • Easternmost (westernmost) feasible location for 
J J satellite j, in degrees west longitude. 


Minimum [E^] over all j (assumed to be zero) 


Maximum [W^] over all j. 


As.. ■ Worst-case minimum required separation between 
satellites i and j, in degrees longitude. 


M = Large positive constant 


J-J 


-M 

if 

A...-0 

As . . 

if 

A, lj > 


Desired location of satellite j, in degrees west 
long i t ud e 

■ (VL - E j ) / 2 (default) 


L . 
J 


Constant inversely proportional to the length of the 
feasible arc of satellite j (used in model 5)* 

- 1 /(W. - E .) . 

J J 
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P. * Population of country with satellite j. 


C. = Constant proportional to the weight assigned to each 
^ satellite j (used in model 9). 

- ( P - ) / I (P.) 

J j J 

Decision Variables: 


x. ■ Orbital location of satellite j (degrees west of E), 
J in degrees west longitude. 


x. + (x.~) * Degrees west (east) of D. that satellite j is 
J ^ located, in degrees we s t ^ longi t ud e . 


x = Dummy satellite location corresponding to the eastern- 
0 most satellite location, in degrees west longitude. 


x * Dummy satellite location corresponding to the 
n+1 westernmost satellite location, in degrees west 
long itud e . 


w. ( e . ) 
J J 


Western (eastern) edge of the arc segment allotted 
satellite j, in degrees west longitude. 


1 if satellite i is located to the west of 
satellite j. 

0 otherwise 



if satellite j is located closer to its west- 
ern boundary than to its eastern boundary. 


0 


otherwise 
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if the western-edge of satellite i is located 
to the west of western-edge of satellite j. 

o th e rwi s e 


n * Minimum gap between adjacent satellites, over all pairs 
of adjacent satellites, in degrees longitude. 


p » Maximum deviation of allotted location from given 
desired satellite Location, over all satellites, in 
degrees longitude. 


q » Maximum deviation of allotted satellite location from 
the closest boundary of satellite's feasible arc, over 
all satellites, in degrees longitude. 


r * Length of the minimum arc segment allotted to a sat- 
ellite, over all satellites, in degrees longitude. 


t - Length of the minimum weighted arc segment allotted to 
a satellite, over all satellites, in degrees longitude. 
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Model 2: Maximize the minimum gap between adjacent 

satellites. 


Maximize 

z = n 



(A.l) 



subject to, 







< x.-Xj) . 

[ (E-W) - A 8 £ j 


> (E-W) 

(A.2) 

Vi, j 

3 i< j 

(-X..X.) - 

[(E-W) -As^ 

'-ii 

> Aa . . 

(A. 3) 

Vi, j 

3 i< j 

( X . - X . ) + 

1 J 

[2(E-W)]x. j - 

11 2 

2 ( E-W ) 

(A. 4) 

Vi, j 

3i<j 

( -X^+x j ) - 

[2(E-W)]x ij - 

n > 

0 

(A. 5) 

Vi, j 

3 i< j 

E . < x . < W 
J - J “ 

j 



( A. 6) 

vj 



(A. 7) Vi,j 3 i< j 
(A. 8) Vj 

n > 0 (A. 9) 

The objective function (A.l) maximizes the minimum gap 
between a pair of adjacent satellites. Constraints (A.2) 
and (A. 3) are dichotomous constraints which ensure that 
minimum satellite separation requirements are enforced. 
Constraints (A. 4) and (A. 5) are also dichotomous 
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constraints and they measure the minimum gap between 
adjacent satellites for all pairs of satellites i and j. 
This minimum gap is denoted by n. Constraint (A. 6) ensures 
that the allotted location for satellites is feasible. The 
integrality of the binary variables is enforced by (A. 7)., 
The non-negativity of the continuous variables is enforced 
by ( A . 8) and ( A. 9 ) • 

For a problem with m satellites, this formulation 
entails 2m(m-l) structural constraints (A. 2), (A. 3), (A. 4) 

and (A. 5), m+1 continuous variables and m(m-l)/2 binary 


variables 
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Model 3: Minimize the total deviation of allotted locations 

from given desired locations, subject to feasible 
arc restrictions [11]. 


Minimize 


z - I ( x. + + x - ) 
i J J 


( A. 10) 


subject to, 


( X £-Xj) + [ ( E- W ) - A 8 ijJ x ij > (E-W) 


(A. 11) Vi, j 3 i< j 




[(E-W) - As ij ]x i . > As. 


( A . 1 2 ) Vi, j 3 i< j 


x . 
J 



x . 


J 


D . 


J 


C A. 13) Vj 


E. < x. < W. (A. 14) Vj 

J J J 

x^ ■ 0 or 1 (A. 15) Vi,j 3i<j 

x^,x^ + ,Xj~ > 0 (A. 16) Vj 

The objective function (A. 10) minimizes the total deviation 
of allotted locations from desired locations. Constraints 
(A. 11) and (A. 12) are dichotomous constraints which ensure 
that minimum satellite separation requirements are 
enforced. The deviation of the allotted location from the 
desired location for each satellite is measured in con- 
straint (A. 13). Constraint (A. 14) ensures that the 
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allotted location for each satellite is feasible* The inte- 
grality of the binary variables is enforced by (A.15)« The 
non-negativity of the continuous variables is enforced by 
( A • 1 6 ) • 

For a problem with m satellites, this formulation 
entails m structural constraints (A. 11), (A. 12) and 

(A.13)> 3m continuous variables and m(m-l)/2 binary vari- 


ables 
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Model 4: Minimize the maximum deviation of allotted location 
from given desired location, subject to feasible 
arc restrictions [20]. 

Minimize z = p (A. 17) 

subject to, 

< x i~x.) + [ < E-W) - A s . . ] x . . > (E-W) (A. 18) Vi,j 3i<j 

( " X i + x j ) ‘ - A8 ij] x ij > As ij (A. 19) Vi , j 3 i< j 

x j - x/ + - D j (A. 20) Vj 

Xj + + Xj" - p < 0 (A. 21) Vj 


E . 
J 


< x j < Wj 


(A. 22) Vj 


Xjj “ 0 or 1 (A. 23) V i , j 3 i< j 

x j,Xj + ,Xj- > 0 (A. 24) Vj 

P > 0 (A. 25) 

The objective function (A. 17) minimizes p, the maximum 
deviation of an allotted location from the corresponding 
desired location. Constraints (A. 18) and (A. 19) are 
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dichotomous constraints which ensure that minimum satellite 
separation requirements are enforced. The deviation of the 
allotted location from the corresponding desired location 
for each satellites is measured in constraint (A. 20)* The 
maximum deviation of an allotted location from its corre- 
sponding desired location is measured by constraint (A. 21). 
Constraint (A. 22) ensures that the allotted location for 
all satellites is feasible. The integrality of the binary 
variables is enforced by (A. 23). The non-negativity of the 
continuous variables is enforced by (A. 24) and (A. 25). 

For a problem with m satellites, this formulation 
entails m(m-t-l) structural constraints (A. 18), (A. 19), 
(A. 20) and (A. 21), 3m+l continuous variables and m(m-l)/2 
binary variables. 
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Model 5: Minimize the total weighted deviation of allotted 

locations from given desired locations [ 2 0 J • 


Minimize 


subject to, 


I L • ( x . + x . ) 

j J J J 


(A. 26) 


( x.- X j) + l(E-W) - As^lx^ > 


(E-W) 


(A. 27) Vi,j 3 i< j 


(-x i + Xj ) - [ (E-W ) -A3 ij ]x ij > As^ (A. 28) Vi,j3i<j 


X . - X . + x . * D . 

J J J J 


(A. 29) Vj 


x £ j - 0 or 1 


(A. 30) Vi,j 3 i< j 


x., x . ,x. > 0 


(A. 31) Vj 


The objective function (A. 26) minimizes the total weighted 
deviation of allotted locations from given desired loca- 
tions. The weight, L ^ , assigned to each satellite j is the 
reciprocal of the length of the feasible arc of satellite 
j. By weighing the objective function in this way we give 
more weight to administrations with smaller arcs and there- 
by attempt to ensure that the solution will be feasible. 
Constraints (A. 27) and (A. 28) are dichotomous constraints 
which ensure that minimum satellite separation requirements 
are enforced. The deviation of the allotted location from 
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the corresponding desired location for each satellite is 
measured in constraint (A. 29). The integrality of the 
binary variables is enforced by (A. 30). The non-negativity 
of the continuous variables is enforced by (A. 31). 

For a problem with m satellites, this formulation 
o 

entails m structural constraints (A. 27), (A. 28) and 

(A. 29), 3m continuous variables and m(m-l)/2 binary vari- 


ables. 
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Model 6: Minimize the maximum deviation of allotted loca- 

tion from the closest boundary of feasible arc. 

Minimize z = q (A. 32) 

subject to, 


( X i-Xj) + [ ( E -W) -As.j] 

x ij > (E-W) 

(A. 33) 

Vi, j 

3i<j 

(_X i + x j ) “ t ( E-W) - As i j ] 

X ij > A °ij 

(A. 34) 

Vi, j 

3 i< j 

(w. - x .) + ( w-e ) x G j - q 

< (W-E) 

(A. 35) 

vj 


U j ‘ E j ) " (W_E)x oj " q 

< 0 

(A. 36) 

vj 


E i < -i < «j 


(A. 37) 

Vj 


x i j = 0 or 1 


( A. 38) 

Vi, j 

3 i< j 

x . = 0 or 1 

oj 


(A. 39) 

Vj 


x j > 0 


( A. 40) 

Vj 


q > 0 


(A. 41) 




The objective function (A. 32) minimizes the maximum devia- 
tion of a satellite's allotted location from the closest 
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boundary of that satellite's feasible arc. Constraints 

(A. 33) and (A. 34) are dichotomous constraints which ensure 
that minimum satellite separation requirements are 
enforced. Constraints (A. 35) and (A. 36) are also dichoto- 
mous or "ei ther-or" constraints. Exactly one of these con- 
straints will be redundant. The binary variable x takes 

o J 

values of 1 or 0, depending on whether satellite j is 

located closer to its western boundary or eastern boundary. 

If satellite j is located closer to its western boundary 

the constraint (A. 36) is redundant and (A. 35) measures the 

deviation from the western boundary. If x .=o then con- 

° J 

straint (A. 36) measures the deviation of satellite j from 
the eastern boundary. The maximum deviation from the clos- 
est boundary is denoted by q. Constraint (A. 37) ensures 
that the allotted location for satellites is feasible. The 
integrality of the binary variables is enforced by (A. 38) 
and (A. 39). The non-negativity of the continuous variables 
is enforced by (A. 40) and (A. 41). 

For a problem with m satellites, this formulation 
entails m(m+l) structural constraints (A 33), (A. 34), 

(A. 35) and (A. 36), m+1 continuous variables and m(m+l)/2 

binary variables. 
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Model 7: Minimize the length of the total arc occupied by 

the satellites to be allotted locations [20]. 


Mi n im i ze 


X n+1 " x 0 


(A. 42) 


subject to, 


( x i" x j) + KE-W) -As.jjx.j > (E-W) 
(-x i+X j) - [(E-W) -As ij Jx i ^ >^®ij 


(A. 43) Vi, j 3 i< j 
(A. 44) Vi, j 3 i< j 
(A. 45) Vj 
(A. 46) Vj 
(A. 47) Vj 


x ij * 0 or 1 (A. 48) Vi , j 3i<j 

x j > 0 (A. 49) Vj 

In this model two dummy variables, x^ and x r + ^, are intro- 
duced. These correspond to the easternmost and westernmost 
allotted satellite locations, respectively. The objective 
function (A. 42) minimizes the total arc occupied by the 
allotted satellites, which is equal to the distance between 
the easternmost and westernmost satellite locations. 
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Constraints (A. 43) and (A. 44) are dichotomous constraints 
which enforce the minimum satellite separation require- 
ments. Constraints (A. 45) and (A. 46) ensure that the dummy 
satellites occupy locations that are separated by at least 
the greatest distance between any two satellites. Con- 
straint (A. 47) ensures that the allotted location for each 
satellite is feasible. The integrality of the binary vari- 
ables is enforced by (A.48)« The non-negativity of the con- 
tinuous variables is enforced by (A. 49)* 

For a problem with m satellites, this formulation 
entails m(m+l) structural constraints (A.43)> (A. 44), 
(A. 45) and (A. 46), m+2 continuous variables and m(m-l)/2 
binary variables. 
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Model 8: Maximize the length of the minimum arc segment 

allotted to any satellite [17]* 

Maximize z = r (A. 50) 

subject to, 


( e . 

l 


w j> 



w i> 


w . 
J 


e . 
J 


e . 
J 


> 


w . 
J 


< "j 


+ [(E-W) 
+ [(E-W) 
r > 0 




> (E-W) 




> 


s. . 

1J 


(A. 51) Vi, j 3 i< j 
(A. 52) V i , j 3 i< j 
(A. 53) Vj 
(A. 54) Vj 
(A. 55) Vj 


= 0 or 1 


(A. 56) Vi, j 3 i< j 


e j > w j 2 0 ( A. 57) Vj 

r > 0 (A. 58) 

The objective function (A. 50) maximizes the length of the 
shortest arc segment allotted to any satellite. Constraints 
(A. 51) and (A. 52) are dichotomous constraints which ensure 
minimum satellite separation requirements by enforcing the 
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required satellite separation between the western edge of 
an allotted arc and the eastern edge of the arc of adjacent 
satellites. By replacing by we permit arcs to 
overlap if the minimum required satellite separation is 
zero. The length of the minimum arc segment allotted to 
any satellite is measured in constraint (A. 53), and is 
denoted by r. Constraints (A. 54) and (A. 55) ensure that 
the arc segment allotted to each satellite lies within its 
feasible arc. The integrality of the binary variables is 
enforced by (A. 56). The non-negativity of the continuous 
variables is enforced by (A. 57) and (A. 58). 

For a problem with m satellites, this formulation 
entails m 4 structural constraints (A. 51), (A. 52) and 
(A. 53), 2tn+ 1 continuous variables and m(m-l)/2 binary vari- 
ables. 


Model 9: Maximize the length of the minimum weighted arc 

segment allotted to any satellite [ 1 7 ] • 
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Maximize 


(A. 59) 


subject to, 


. Wj ) + l(E-W) - Sijlx.j > 


(E-W) (A. 60) Vi, j 3 i< j 


(e - w.) + [(E-W) - i.jjx.j > i-j (A. 61) Vi,j 3i<j 


w. - e . - c . t > 0 
J J J ~ 


(A. 62) Vj 


e j * E j 


(A. 63) Vj 


"j 1 W J 


(A. 64) Vj 


= 0 or 1 


(A. 65) Vi, j 3 i< j 


e., w . > 0 


(A. 66) Vj 


t > 0 


(A. 67) 


The objective function (A. 59) maximizes the length of the 
shortest weighted arc segment allotted to any satellite. 
Constraints (A. 60) and (A. 61) are dichotomous constraints 
which ensure minimum satellite separation requirements by 
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enforcing the required separation between the western edge 
of an allotted arc and the eastern edge of the arc of adja- 
cent satellites. By replacing As. .by b- we permit arcs 

^ J ^ J 

to overlap if the minimum required satellite separation is 
zero. The length of the minimum weighted arc segment 
allotted to any satellite is measured in constraint (A. 62), 
and is denoted by t. Constraints (A. 63) and (A. 64) ensure 
that the arc segment allotted to each satellite lies within 
its feasible arc. The integrality of the binary variables 
is enforced by (A.65)« The non-negativity of the continuous 


variables is enforced by (A. 66) and (A. 67) 


For a problem with m satellites, this formulation 
entails m^ structural constraints (A. 60), (A. 61) and 
(A. 62), 2m+l continuous variables and m(m-l)/2 binary vari- 
ables. 


Appendix B 

COMPUTER RESULTS OF MIP FORMULATIONS 


84 



85 


Table 11 

COM PU TER R ESULTS OF M1P FORMULATIONS ( E. E U RO PE) 


FIRST FEASIBLE 
SOLUTION 

M on F T . 

BEST SOLUTION 

FINAL SOLN. 

nuuL L — ■ — 




Obj Fn . Time 

Obj Fn. 

Time 

Time 

Value Min 

Value 

Min. 

Min. 


1 

4.54 

.31 

9.33 

3.45 

CM 

r^ 

O 

• 

O 

CM 

2 

7.45 

.71 

11.27 

1.66 

20 . 20 2 

3 

77.06 

.15 

49.87 

10.21 

20. 12 2 

4 

12.32 

.19 

8.48 

18.95 

20 .06 2 

5 

69.75 

.13 

58.21 

3.53 

20. 03 2 

6 

14.22 

.11 

2.71 

3.58 

20 .09 2 

7 

32.81 

.23 

25.48 

7.86 

20. 21 2 

00 

5.51 

.19 

8.75 

14.49 

20. 22 2 

9 

6.43 

.17 

10.06 

1 .46 

20. 19 2 


NOTE: 1 - Proven Optimum solution 

2 - Terminating at the maximum allotted CPU time. 
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Table 12 



COMPUTER 

RESULTS OF 

MIP FORMULATIONS (W. 

EUROPE) 

MODEL 

FIRST FEASIBLE 
SOLUTION 

BEST 

SOLUTION 

FINAL SOLN. 


Ob j Fn 
Value 

. Time 

Min • 

Ob j Fn 
Value 

• Time 

Min* 

Time 

Min. 

i 

4.77 

.29 

8.38 

10.38 

20 . II 2 

2 

8.22 

.55 

10.36 

1 . 68 

20 • 1 8 2 

3 

61.92 

.09 

29.67 

0.69 

5.22 1 

4 

7.25 

• 1 7 

6 .62 

0.56 

2 . 5 9 1 

5 

46.12 

.10 

26.69 

3.45 

7 • 1 5 1 

6 

12.33 

.17 

7.51 

18.67 

20. II 2 

7 

33.13 

.26 

23.71 

19.35 

20 . 02 2 

8 

5.61 

.10 

6.00 

0.21 

0 . 2 8 1 

9 

7 .79 

.21 

10.14 

19.73 

20 . 06 2 


NOTE: 1 - Proven Optimum solution 


2 - Terminating at the maximum allotted CPU time 
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Table 13 



COMPUTER 

RESULTS OF 

MIP FORMULATIONS (S. 

AMERICA) 






MODEL 

FIRST FEASIBLE 
SOLUTION 

BEST SOLUTION 

FINAL SOLN. 

Ob j Fn 
Value 

. Time 

Min. 

Obj Fn . 
Va lue 

Time 

Min. 

Time 

Min. 

1 

7.77 

.41 

10.55 

11.83 

20. 22 2 

2 

11.27 

1 .05 

13.20 

15.25 

20. 15 2 

3 

33.55 

.13 

30.44 

1.44 

20. 02 2 

4 

6.96 

.19 

6.60 

0.21 

20. 02 2 

5 

38.02 

.10 

27.70 

1.30 

20. 04 2 

6 

41.64 

.09 

0.86 

17.07 

20. 05 2 

7 

28.66 

.35 

22.50 

6.27 

20. 08 2 

8 

7 .68 

.37 

9.85 

16.41 

20. 24 2 

9 

9.67 

.15 

11.99 

2.30 

20. 01 2 

NOTE : 

1 - Proven Optimum 

so lu t ion 




2 - Terminating at the maximum allotted CPU time. 
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Table 14 

COMPUTER RESULTS OF MIP FORMULATIONS (N. AFRICA) 







— — — 



FIRST 

FEASIBLE 

BEST 

SOLUTION 

FINAL SOLN. 


SOLUTION 




MODEL— 



— 

— 

— 


Ob j Fn 

. Time 

Ob j Fn 

• Time 

Time 


Value 

Min 

Value 

Min. 

Min . 

i 

11.10 

.15 

16.90 

2.11 

20. 13 2 

2 

12.67 

.30 

18.44 

5.36 

CM 

O 

CM 

• 

O 

CM 

3 

8.65 

.01 

8.65 

.01 

0.06 1 

4 

1 • 64 

.01 

1 . 64 

.01 

0 . 0 7 1 

5 

6.85 

.01 

6.85 

.01 

0 . 05 1 

6 

2 . 99 

.01 

0.00 

.77 

0.81 1 

7 

17.48 

.10 

16.27 

15.47 

20. 04 2 

8 

11.35 

.15 

14.85 

8.36 

CM 

O' 

^4 

m 

o 

CM 

9 

12 . 37 

.07 

14.42 

4.27 

CM 

O 

• 

o 

CM 


NOTE: 1 - Proven Optimum solution 

2 - Terminating at the maximum allotted CPU time. 
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Table 15 

COMPUTER RESULTS OF MIP FO RMULATIONS (S.E. ASIA) 


MODEL 

FIRST FEASIBLE 
SOLUTION 

Obj Fn . Time 

Value Min 

BEST 

Obj Fn 
Va lue 

SOLUTION 

• Time 

Min. 

FINAL SOLN. 

Time 

Miri. 

1 

9.69 

.13 

14.69 

17.14 

CNI 

Os 

O 

• 

o 

CM 

2 

14.00 

.26 

16.44 

.75 

CM 

O 

CM 

• 

O 

CM 

3 

24.11 

.02 

23.05 

.54 

1 .60 1 

4 

5.88 

.05 

5.19 

. 1 1 

1 . 33 1 

5 

26.59 

.03 

18.31 

.14 

0. 91 1 

6 

12.25 

.05 

0.74 

6.05 

20 . 04 2 

7 

26.16 

.08 

18.99 

18.88 

CM 

o 

. 

o 

CM 

8 

9.06 

. 1 1 

13.18 

6.93 

20. li 2 

9 

8.67 

.03 

8.67 

.03 

0 . 09 1 

NOTE: 1 

- Proven 

Op t imum 

solution 




2 - Terminating at the maximum allotted CPU time. 




COMPARISON OF SOLUTION QUALITY OF MIP MODELS 


91 


Table 16 

COMPA RISON OF SOLUTION Q UALI TY ( E. EU RO PE) 


BEST 

SOLN. OBJECTIVE FUNCTION VALUE IN MODEL 

TO 


MODEL 

1 

2 

3 

4 

5 

6 

7 

1 

9.3 

9.3 

3 76.8 

62.0 

357 . 7 

55.7 

124.0 

2 

6.7 

11.3 

412.8 

62.0 

404 . 9 

50.4 

124.0 

3 

0.0 

0.6 

49.9 

10.7 

51.2 

61 .0 

28.4 

4 

0.0 

0.5 

71.6 

8.5 

70.2 

54.3 

27.3 

5 

0.0 

0.3 

60.9 

14.5 

58.2 

59.5 

32.0 

6 

0.0 

1.3 

593.0 

61.0 

563.1 

2 . 7 

124.0 

7 

0.0 

0.4 

270.6 

38.4 

282.3 

47.9 

25.5 




c 
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Table 17 

C OMPARI S ON OF SOLUTION QUALITY (W. EUROPE) 


BEST 

SOLN. 

TO 

MODEL 

1 

OBJECTIVE FUNCTION VALUE IN MODEL 


2 

3 

4 

5 

6 

7 

i 

8.4 

9.3 

349.5 

59.0 

365.4 

55.5 

114.0 

2 

8.0 

10.4 

366.6 

59.0 

358.3 

50.5 

114.0 

3 

0.0 

1.3 

29*7 

11.8 

26.6 

61.4 

34 . 8 

4 

0.0 

1.8 

53.4 

6.6 

96.7 

61.4 

35.3 

5 

0.0 

1.3 

29.7 

11.8 

26.6 

61.4 

34 . 8 

6 

0.0 

1.3 

575.9 

62.0 

568.0 

7.5 

120.0 

7 

0.0 

0.9 

230.6 

35.3 

256.1 

53.1 

23 . 7 


93 


Table 18 

C OMPARI SON OF SOLU TI ON QUAL ITY (S. AMERICA ) 


BEST 

SOLN. 

TO 

MODEL 

1 

OBJECTIVE FUNCTION VALUE IN MODEL 

2 3 4 5 6 

7 

1 

10.6 

10.6 

496.0 

68.0 

407.3 

61.4 

152.0 

2 

8.1 

13.2 

530.4 

68.0 

436.5 

55.6 

160.0 

3 

0.0 

0 • 6 

30.4 

1 1.4 

28.4 

68.0 

35.0 

4 

0.0 

0.2 

55.3 

6.6 

46.3 

66.8 

22.5 

5 

0.0 

0.3 

30.4 

10.0 

27.7 

67.3 

33 . 8 

6 

0.0 

0.5 

7 74.7 

68.0 

647.2 

0.9 

160. 0 

7 

0.0 

0.0 

556.4 

57.6 

468 . 7 

29.0 

22.5 
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Table 19 

COMPARISON OF SOLUTION QUALITY ( N ♦ A FRICA) 


BEST 

SOLN. 

TO 

MODEL 

1 

OBJECTIVE FUNCTION VALUE IN MODEL 
2 3 4 5 6 

7 

1 

16.9 

1 7.4 

400 . 7 

64.0 

323.5 

57.8 

166.0 

2 

18.4 

13.4 

431.1 

63.0 

350.8 

45.8 

166.0 

3 

0 . 0 

2.0 

8.7 

3.3 

7 . 1 

64.3 

43.0 

4 

0 . 0 

0.7 

13.0 

1 . 6 

10.4 

63.4 

44.3 

5 

0.0 

2.0 

8.7 

3 . 3 

6.9 

64.0 

46.3 

6 

2 . 9 

6 0 

616.0 

65.0 

500.0 

0.0 

166.0 

7 

0 . 0 

0.9 

363.4 

60.0 

295.3 

48.8 

16.3 
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Table 20 

COMPARISON OF S OLUTIO N QUALITY (S.E. ASIA) 


BEST 

SOLN. 

TO 

MODEL 


OBJECTIVE 

FUNCTION 

VALUE IN 

MODEL 


1 

2 

3 

4 

5 

6 

7 

1 

14.7 

15-3 

386.9 

64.0 

316.3 

47.1 

148.0 

2 

12.2 

16.4 

391.1 

64.0 

314.5 

49.8 

148.0 

3 

0.0 

1 .4 

23.0 

7.1 

18.4 

66 . 7 

27.6 

4 

0.0 

0.4 

40.0 

5-2 

34 . 6 

62.8 

32.2 

5 

0.0 

1 .4 

23.1 

9.8 

18.3 

66 . 8 

27.8 

6 

0.4 

1.4 

588.7 

67.0 

498.9 

0.7 

144.0 

7 

0.0 

0.0 

319.9 

47 . 8 

279.8 

37.9 

19.0 


Appendix D 


ORDERING OF SATELLITE ALLOTMENTS IN MIP MODELS 
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Table 21 

CHANGE IN THE ORDERING OF SATELLITES (E. EUROPE) 


OBJ 

1(F) 

1 

8 

10 

6 

4 

12 

9 

1 1 

7 

5 

3 

2 

OBJ 

1(B) 

1 1 

7 

12 

3 

10 

8 

1 

2 

6 

4 

9 

5 

OBJ 

2(F) 

1 

10 

1 1 

12 

6 

9 

8 

7 

4 

3 

5 

2 

OBJ 

2(B) 

8 

3 

11 

6 

10 

2 

9 

1 

12 

7 

4 

5 

OBJ 

3 ( F) 

2 

3 

1 

4 

7 

10 

9 

1 1 

12 

8 

6 

5 

OBJ 

3(B) 

1 1 

12 

8 

9 

1 

5 

7 

10 

4 

6 

2 

3 

OBJ 

4 ( F) 

2 

3 

1 

4 

9 

1 1 

12 

8 

10 

7 

5 

6 

OBJ 

4(B) 

12 

8 

6 

5 

1 1 

10 

2 

9 

1 

7 

4 

3 

OBJ 

5 ( F) 

2 

3 

1 

5 

7 

10 

9 

8 

6 

4 

1 1 

12 

OBJ 

5(B) 

9 

12 

11 

6 

4 

1 

7 

5 

1 0 

3 

8 

2 

OBJ 

6(F) 

2 

5 

4 

3 

11 

7 

9 

6 

1 

10 

12 

8 

OBJ 

6(B) 

12 

5 

7 

2 

6 

1 

10 

1 1 

9 

8 

3 

4 

OBJ 

7(F) 

3 

9 

10 

12 

4 

7 

1 1 

8 

6 

5 

1 

2 

OBJ 

7(B) 

1 

8 

9 

4 

6 

2 

1 1 

10 

7 

5 

12 

3 

OBJ 

8(F) 

2 

5 

1 

4 

6 

1 1 

10 

8 

9 

12 

7 

3 

OBJ 

8(B) 

2 

4 

11 

7 

1 

5 

9 

10 

3 

6 

8 

12 

OBJ 

9(F) 

2 

5 

4 

1 

10 

3 

6 

7 

12 

9 

1 1 

8 

OBJ 

9(B) 

3 

7 

4 

9 

2 

1 

12 

5 

10 

8 

6 

1 1 
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Table 2 2 

CHANGE IN THE ORDERING OF SATELLITES (W. E U RO P E) 


OBJ 

1(F) 

3 

2 

8 

9 

10 

11 

12 

7 

4 

5 

6 

1 

OBJ 

KB) 

1 

9 

12 

4 

11 

6 

8 

2 

7 

10 

3 

5 

OBJ 

2(F) 

5 

3 

2 

9 

1 1 

12 

8 

10 

6 

7 

4 

1 

OBJ 

2(B) 

1 

9 

7 

4 

11 

12 

8 

2 

10 

6 

3 

5 

OBJ 

3(F) 

2 

3 

1 

9 

4 

6 

7 

5 

10 

1 1 

12 

8 

OBJ 

3(B) 

1 1 

12 

7 

10 

9 

8 

6 

5 

3 

1 

4 

2 

OBJ 

4 ( F) 

2 

3 

1 

4 

5 

6 

10 

7 

8 

12 

9 

1 1 

OBJ 

4(B) 

1 1 

9 

12 

8 

7 

10 

4 

6 

5 

3 

1 

2 

OBJ 

5( F) 

2 

1 

3 

5 

6 

7 

8 

4 

10 

11 

12 

9 

OBJ 

5(B) 

1 1 

12 

7 

10 

9 

8 

6 

5 

3 

1 

4 

2 

OBJ 

6(F) 

5 

1 

4 

7 

6 

1 1 

12 

9 

2 

3 

8 

1 0 

OBJ 

6(B) 

7 

1 

4 

3 

2 

9 

1 1 

10 

8 

12 

6 

5 

OBJ 

7(F) 

3 

5 

6 

10 

12 

7 

8 

9 

1 1 

4 

2 

1 

OBJ 

7(B) 

4 

8 

12 

9 

1 

3 

11 

5 

2 

6 

10 

7 

OBJ 

8(F) 

4 

1 

6 

7 

2 

5 

1 2 

3 

10 

8 

1 1 

9 

OBJ 

8(B) 

7 

4 

1 1 

1 

2 

9 

6 

3 

5 

12 

8 

1 0 

OBJ 

9(F) 

4 

5 

3 

7 

2 

8 

12 

11 

9 

10 

6 

1 

OBJ 

9(B) 

5 

4 

1 1 

6 

2 

7 

12 

3 

1 

8 

9 

1 0 
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Table 23 

CHANGE IN THE ORDERING OF SATELLITES (S» AM E RI C A) 


OBJ 

1 ( F) 

2 

13 

9 

5 

4 

12 

7 

3 

10 

11 

1 

8 

6 

OBJ 

1(B) 

9 

7 

4 

13 

1 

5 

11 

12 

8 

6 

3 

2 

10 

OBJ 

2(F) 

2 

3 

13 

4 

12 

5 

11 

8 

1 

9 

7 

10 

6 

OBJ 

2(B) 

6 

9 

1 

8 

4 

7 

13 

2 

3 

5 

12 

11 . 

10 

OBJ 

3(F) 

3 

13 

10 

1 1 

8 

7 

2 

12 

1 

5 

4 

9 

6 

OBJ 

3(B) 

6 

9 

1 

5 

4 

12 

2 

7 

8 

11 

10 

13 

3 

OBJ 

4 ( F) 

3 

13 

10 

8 

7 

11 

2 

12 

1 

5 

4 

9 

6 

OBJ 

4(B) 

6 

9 

4 

5 

1 

12 

2 

7 

1 1 

13 

8 

10 

3 

OBJ 

5(F) 

3 

1 1 

13 

10 

8 

7 

2 

1 

5 

4 

12 

9 

6 

OBJ 

5(B) 

6 

9 

1 

12 

4 

5 

2 

7 

8 

1 1 

10 

13 

3 

OBJ 

6 ( F) 

1 0 

11 

2 

8 

3 

13 

9 

6 

12 

5 

4 

7 

1 

OBJ 

6(B) 

6 

2 

5 

7 

1 

13 

12 

4 

9 

3 

8 

1 1 

10 

OBJ 

7(F) 

3 

12 

1 1 

13 

10 

9 

8 

7. 

4 

6 

2 

5 

1 

OBJ 

7(B) 

13 

1 

9 

12 

8 

6 

1 1 

7 

2 

5 

10 

4 

3 

OBJ 

8(F) 

2 

5 

4 

6 

11 

7 

9 

10 

12 

13 

3 

8 

1 

OBJ 

8(B) 

10 

2 

11 

12 

8 

6 

13 

3 

1 

5 

4 

7 

9 

OBJ 

9(F) 

7 

1 1 

10 

8 

3 

6 

1 

5 

13 

12 

2 

4 

9 

OBJ 

9(B) 

7 

8 

10 

11 

3 

13 

4 

12 

1 

5 

2 

9 

6 
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Table 24 



CHANGE 

IN 

THE 

ORDERING 

; of 

SATELLITES 

(N . 

AFRICA) 

OBJ 

1(F) 

10 

5 

1 

4 

6 

3 

8 

7 9 

OBJ 

1(B) 

5 

2 

4 

1 

3 

7 

10 

8 9 

OBJ 

2(F) 

10 

5 

4 

1 

6 

7 

8 

9 3 

OBJ 

2(B) 

5 

3 

1 

9 

8 

4 

10 

2 7 

OBJ 

3(F) 

6 

7 

8 

1 

9 

2 

10 

3 4 

OBJ 

3(B) 

5 

4 

3 

10 

2 

9 

1 

8 7 

OBJ 

4(F) 

6 

7 

8 

1 

9 

2 

10 

3 4 

OBJ 

4(B) 

5 

4 

3 

10 

9 

2 

1 

8 7 

OBJ 

5 ( F) 

6 

7 

8 

1 

9 

2 

10 

3 4 

OBJ 

5(B) 

5 

4 

3 

10 

2 

9 

1 

8 7 

OBJ 

6 ( F) 

6 

7 

8 

1 

9 

2 

10 

3 5 

OBJ 

6 ( B) 

5 

3 

9 

8 

7 

4 

10 

2 1 

OBJ 

7(F) 

2 

4 

7 

8 

9 

5 

10 

6 3 

OBJ 

7(B) 

1 

5 

7 

1 0 

8 

3 

9 

6 4 

OBJ 

8(F) 

1 

3 

5 

6 

10 

7 

4 

8 9 

OBJ 

8(B) 

7 

9 

8 

1 0 

6 

3 

1 

4 2 

OBJ 

9 ( F) 

7 

2 

8 

1 0 

6 

9 

4 

1 5 

OBJ 

9(B) 

7 

2 

6 

10 

1 

3 

8 

9 4 


2 

6 

2 

6 

5 

6 

5 

6 

5 

6 

4 
6 

1 

2 

2 

5 

3 


5 
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Table 25 



CHANGE 

IN 

THE 

ORDERING 

OF 

SATELLITES 

( S .E 

. ASIA) 



OBJ 

1(F) 

1 

9 

6 

4 

7 

8 

10 

5 

3 

2 

OBJ 

1(B) 

10 

4 

7 

6 

3 

8 

5 

1 

9 

2 

OBJ 

2(F) 

1 

6 

4 

8 

10 

5 

3 

9 

7 

2 

OBJ 

2(B) 

10 

5 

2 

6 

7 

8 

3 

9 

4 

1 

OBJ 

3(F) 

3 

1 

2 

4 

7 

8 

5 

9 

6 

10 

OBJ 

3(B) 

10 

6 

9 

5 

8 

7 

4 

1 

2 

3 

OBJ 

4(F) 

3. 

1 

7 

2 

4 

5 

8 

6 

9 

10 

OBJ 

4(B) 

10 

9 

4 

6 

5 

8 

7 

3 

2 

1 

OBJ 

5(F) 

1 

3 

2 

4 

7 

8 

5 

10 

9 

6 

OBJ 

5(B) 

9 

10 

6 

5 

8 

7 

4 

1 

2 

3 

OBJ 

6 ( F) 

2 

5 

4 

7 

6 

9 

3 

10 

8 

1 

OBJ 

6(B) 

9 

4 

7 

1 

8 

10 

3 

6 

5 

2 

OBJ 

7(F) 

6 

1 

4 

9 

8 

10 

7 

5 

2 

3 

OBJ 

7(B) 

9 

6 

3 

8 

5 

2 

7 

10 

1 

4 

OBJ 

8(F) 

7 

9 

1 

4 

10 

8 

5 

3 

6 

2 

OBJ 

8(B) 

2 

3 

6 

7 

8 

5 

9 

1 

4 

10 

OBJ 

9(F) 

1 

2 

5 

3 

4 

9 

8 

6 

7 

10 

OBJ 

9(B) 

1 

2 

5 

3 

4 

8 

9 

6 

7 

10 
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